Transport in Coupled Quantum Dots: Kondo Effect Versus Anti-Ferromagnetic 
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The interplay between the Kondo effect and the inter-dot magnetic interaction in a coupled-dot 
system is studied. An exact result for the transport properties at zero temperature is obtained by 
diagonalizing a cluster, composed by the double-dot and its vicinity, which is connected to leads. It 
is shown that the system goes continuously from the Kondo regime to an anti-ferromagnetic state 
as the inter-dot interaction is increased. The conductance, the charge at the dots and the spin-spin 
correlation are obtained as a function of the gate potential. 
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In the last years electron transport through a quan- 
tum dot (QD) has been the subject of many experimen- 
tal and theoretical investigations. The possibility of 
measuring the Kondo effect, the basis of heavy- fermion 
physics, in a system with a QD has been proposed by 
several theoretical studies pLE}. Its detection is a dif- 
ficult task since it depends on several different energy 
scales and their relative sizes, such as the coupling con- 
stant between the QD and leads, the Kondo tempera- 
ture Tr and the energy spacing between the dot levels, 
which implies small size dots. Such a very small QD was 
recently obtained and the Kondo effect was measured 
|^| . This confirms the possibility of studying many of the 
properties of strongly correlated metals and insulators in 
artificially constructed mesoscopic structures. The ad- 
vantage of these systems is that their parameters can be 
continuously changed by modifying the applied external 
potentials so that different regimes can be studied. 

The physics associated to heavy-fermion compounds in 
the vicinity of what is called the quantum phase transi- 
tion is signed by the competition between the spin corre- 
lation among the magnetic atoms and between the atoms 
and the conduction electrons (t| . The latter gives rise to 
the Kondo effect and the former tends to create a ferro 
or anti-ferromagnetic ground state destroying the Kondo 
regime. The study of the two impurity Kondo Hamilto- 
nian was proposed to clarify this important problem, as it 
includes in the model the two interactions responsible for 
this competing behavior The competition gives rise 
to several phenomena such as magnetic impurity corre- 
lation, one and two-stage Kondo effect and, according to 
the particle-hole symmetry, Fermi and non-Fermi liquid 
behavior ||]. 

There have been extensive experimental transport 
studies in the double-dot system ju| . They were mainly 
designed to analyze the double-dot molecule. From the 
theoretical view point transport through a double-dot 
has received a considerable attention mainly restricted to 
the study of high temperature, Coulomb blockade phe- 
nomenon uM, 



At very low temperature, the study of a completely 
different physics, created by the interplay of the Kondo 
effect and the inter-dot anti-ferromagnetic correlation, is 
now feasible due to the recent possibility of constructing 
very small dots. An interesting analysis of these phenom- 
ena in a two-QD system, based on qualitative arguments 
and slave-boson mean-field theory, has recently appeared 
||. The behavior of the conductance was studied for a 
system where the double occupancy at the dots was elim- 
inated from the Hilbert space. Although this method be- 
comes exact when the number of spin degrees of freedom 
is infinite, some caution is necessary when applying it to 
spin 1/2 systems. 

The purpose of this letter is to present a numerically 
exact calculation of the transport properties of a double- 
dot system to investigate the competition between the 
Kondo effect and the anti-ferromagnetic correlation. The 
parameters are chosen to reflect real experimental condi- 
tions and, in particular, the intra-dot Coulomb repulsion 
is taken to be finite. The conductance is obtained for sev- 
eral inter-dot coupling constants. The density of states, 
the charge inside the dots and the various spin-spin cor- 
relation functions are calculated in order to characterize 
the state of the system. 

The very low temperature properties (T <C Tk ) of the 
system are obtained by using the Lanczos method Jl^] 
to calculate the ground state of a small cluster contain- 
ing the two dots, which is embedded into the leads. The 
conductance calculated according to this procedure tends 
very rapidly to its exact result as the size of the cluster 
is increased. From this view point our calculation is nu- 
merically exact since we were able to reach convergence, 
within 1% error, for a small cluster. For T > Tk we 
use a self-consistent solution for the equation of motion 
for the one-particle Green function |Dj ] . This approxima- 
tion eliminates all the low lying excitations in the vicinity 
of the Fermi level and, as a consequence, the Kondo ef- 
fect. However, it provides an adequate description of the 
Coulomb blockade high temperature regime ||. 

The system is represented by an Anderson two- 
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impurity first-neighbor tight-binding Hamiltonian, 

H = H c + tY / c+c ja (1) 

where t is the nearest-neighbor hopping in the leads. The 
Hamiltonian H c for the cluster of 2M + 2 sites contain- 
ing the two dots denoted by a and /3 and the other sites 
numbered from 1 to M and 1 to M can be written as 

H c = V ^ "™ + ^ E n ran r g + t" ^(c+^C^ + C.C.) 

+t' X)(4 ct c Ict + c+ a c la + c.c.) + t ]T c+c ja (2) 

er 

where U and V represent, respectively, the electronic re- 
pulsion and the gate potential at the dots, t" is the inter- 
dot interaction, if the coupling between the dots and the 
leads, and the sub- indices i and j, in the last term, run 
over the sites of the cluster other than the dot sites. We 
restrict our analysis to the case of two identical dots. 
An inter-dot magnetic exchange term J is not explicitly 
included in the Hamiltonian since it is not an indepen- 
dent parameter but a function of the inter-dot tunneling 

(j~t»yu)§. 

To describe the very low temperature physics of the 
system, we calculate the one particle Green functions 
G a p at the dots. They are imposed to satisfy a Dyson 
equation G — g + gTG where g is the cluster Green func- 
tion matrix and T is the matrix of the coupling Hamil- 
tonian between the cluster and the rest of the system. 
The undressed Green function g is calculated using the 
cluster ground state obtained by the Lanczos method. 
The Dyson equation proposed is equivalent to the chain 
approximation in a kinetic energy diagrammatic expan- 
sion for the Hubbard Hamiltonian [14]. To be consis- 
tent, the charge of the dressed and undressed cluster 
is imposed to be the same. We calculate g as a com- 
bination of the Green functions of n and n + 1 elec- 
trons with weights (1 — p) and p, g = (1 — p)g n + pg n +i 
M|. In this case the charge of the undressed cluster is 
q c = (1— p)n+p(n+l). The charge of the cluster linked to 
the leads can be expressed as Q c = 2 jy ^ Im.Gu(tu)dtu 
where i runs through all the cluster sites and the factor 
2 includes the electronic spin. This equation together 
with the imposed condition q c = Q c constitute a sys- 
tem of two equations, which requires a self-consistent so- 
lution to obtain p and n as a function of the parame- 
ters of the system. Due to the symmetry of the prob- 
lem, it is useful to represent the Hamiltonian using the 
two-dot bonding and anti-bonding states. In this case, 
according to the Keldysh formalism Jig] , the conduc- 
tance can be written as a = e 2 / h[t' 2 p(e j)] 2 |Gj — G^\ 2 
where G^ = 1/(lo — e ± — S c (w) — T,^(w)) are the one 
particle Green functions at the dots represented in the 
bonding and anti-bonding wave functions, with energies 



e ± = V±t' . S c (w) and Y£(w) are the self-energies corre- 
sponding to the leads and to the many-body interaction, 
respectively. 

The expression for the current is obtained by con- 
sidering explicitly that the system is a Fermi liq- 
uid lim w „»o I m ^J(w) oc lu 2 , which implies that 
there are no dissipative processes taking place at 
the Fermi level. It is important to point out that 
our diagonalization self-consistent procedure satisfies, 
without any imposition, the Luttinger-Ward identity 
Im fJ^dY^iuj) / du)G±{uj)dLU = @. This ensures 
the fulfillment of the Fridel Sum Rule and of the Fermi 
liquid properties. 

We have as well studied this problem through the 
equation of motion for the G + and G~ using a self- 
consistent decoupling procedure, proposed by Hubbard 
|l3| . Since it eliminates all the low lying excitations 
involved in the Kondo effect and properly describes 
the high energy states we suppose it to be an ade- 
quate solution for high temperatures, T > T^. It has 
been used with success to study the Coulomb block- 
ade behavior of the current going through a highly con- 
fined region, as in a QD g). Within this approxi- 
mation the many-body self-energy can be expressed as 
S±( w ) = n±[/(o; - e±)/(w - ir\ ~ e ± - (1 - n ± )J7). It 
requires the self-consistent calculation of the charge at 
the bonding and anti-bonding state ri^ . As shown below 
the two treatments are equivalent in the limit of strong 
anti-ferromagnetic inter-dot coupling regime, where there 
is no Kondo ground state. 

In order to provide a complete characterization of 
the ground state we have compared the inter-dot spin- 
correlation {SdaSdp) with the spin-correlation of one dot 
and its nearest-neighbor conduction electrons {SdSdp). 
For the sake of simplicity, we have calculated these ob- 
jects on the cluster ground state. 

The transport properties are studied for weak and 
strong inter-dot coupling. We solve an eight-atom clus- 
ter with the parameter T/U = t' 2 /WU = 0.08, where W 
is the leads bandwidth. This value corresponds approxi- 
mately to the device where the Kondo effect has recently 
been measured and is compatible with first princi- 
ple calculated parameters The current, the dots 
charge and the spin-spin correlation for weak inter-dot 
coupling (t" /U — 0.16), as a function of V, are presented 
at the left part of Fig. |l|. For V > 0.2 the dot levels are 
above the Fermi level and there is no current. As V is 
reduced, the two dots become partially charged because 
the anti-bonding state enters into resonance. When the 
total charge at the dots is N ~ 1 there is a peak in the 
conductance which reaches the value of e 2 /h. As soon 
as the system enters into the fluctuating valence regime, 
corresponding to e~ below the Fermi energy, the Kondo 
peak appears creating a new channel for the electrons 
to flow. These features are shown in Figs, la and lb. 
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The spin correlation {S c iSdp), continuous line in Fig. lc, 
increases. As V is further augmented extra charge en- 
ters into the dots and the conductance diminishes. As 
e~ goes far beyond the Fermi energy the Kondo temper- 
ature reduces exponentially so that t" 2 /(UTk) increases, 
partially quenching the Kondo effect. The conductance 
has a minimum for V/U = —0.5, when the system has 
electron-hole symmetry and N ~ 2. However, the current 
is large since the system is still within the Kondo regime, 
as reflected by the relationship (SdaSd/i) *C (SdSdp), 
shown in Fig.lc. The large conductance in the inter- 
peak region is due to the Abrikosov-Suhl resonance cre- 
ated by the Kondo effect. This channel is absent when 
T > Tk- The Hubbard solution, also shown in Fig. la, 
can be taken to be a high temperature description of the 
same problem. It preserves the structure of the Coulomb 
blockade, with no conductance between the peaks, which 
are splitted into two due to the inter-dot interaction. 
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FIG. 1. Conductance, charge and spin correlation with 
V/U = 0.08, for weak (t"/U = 0.16) and strong (t"/U = 2.0) 
coupling. Conductance in units of e 2 /h: weak (a) and strong 
(d) coupling; exact-cluster (continuous) and Hubbard approx- 
imation (dotted). The inset in (d) corresponds to the interme- 
diate coupling {t" /U = 0.4) in the exact-cluster approxima- 
tion. Charge: weak (b) and strong (e) coupling; exact-cluster 
results for charge per spin in each dot(continuous) and in 
the +— states (dotted); Hubbard approximation results for 
charge per spin in each dot (dashed). Spin correlation: 
weak (c) and strong (f) coupling; dot-conduction-electron, 
(SciSdp), (continuous) and dot-dot, (SdaSdp), (dotted). 



It is important to notice the different way the charge 
enters into the dots in the Coulomb blockade and in the 
Kondo regimes. In the former, the charge presents a 
plateau due to electronic repulsion while in the latter, 
the Kondo peak going through the Fermi level permits 
a continuous entrance of the charge. This behavior is 
shown in Fig. lb. 

Finally in the case of the strong inter-dot interaction 
(t"/U = 2), shown at the right part of Fig. 1, the Kondo 
regime is almost quenched by the anti-ferromagnetic cou- 
pling. The Kondo spin-spin correlation, although small, 
has its maximum in the near vicinity of the conductance 
peaks, where the Kondo temperature is high, while in 
the whole intermediate region anti-ferromagnetism con- 
trols the conduction. In this case, the conductance re- 
sults obtained by the two approaches, shown in Fig. Id, 
are almost identical, reflecting the fact that the system 
is outside the Kondo regime for most values of the gate 
potential. In the region between the peaks the Hubbard 
approximation presents a very small contribution to the 
conductance, due to the inter-dot interaction splitting, 
which is absent in the exact-cluster solution. 

When the dot-dot interaction has an intermediate 
value (t"/U — 0.4) the conductance in the region be- 
tween the two peaks results to be slowly dependent upon 
the gate potential, as shown in the inset of Fig. Id. This 
behavior has been characterized as a two-plateau struc- 
ture in the conductance as a function of the gate voltage 
by a mean-field slave-boson approximation || , which gets 
a strictly zero conductance in the vicinity of the gate po- 
tential value for which the system has electron-hole sym- 
metry. However we do not get neither a strict plateau 
structure nor a complete quenching of the Kondo effect. 
In the inter-peak region, for finite Coulomb repulsion U 
the current is not zero, although it is small when the dots 
are strongly coupled. The dependence of the conductance 
for the electron-hole symmetric situation as a function of 
the inter-dot coupling is shown in Fig. 2a. It is clearly 
seen that the interplay between the Kondo effect and 
the anti-ferromagnetic correlation does not produce any 
abrupt transition. Our results show a crossover behavior 
as the inter-dot interaction is varied. 




FIG. 2. Conductance in units of e 2 /h (a) as a func- 
tion of the inter-dot coupling for V/U = —0.5, (b) as a 
function of gate potential for cluster of 4 sites(dashed), 8 
sites (continuous), 12 sites (crosses) 
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In order to study the convergence of our results with 
the size of the cluster we present in Fig. 2b the weak cou- 
pling limit conductance for three cluster sizes. A compar- 
ison of the conductance for eight and twelve-site clusters, 
which differs at most by 1%, allows us to conclude that 
our results are numerically exact. 
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FIG. 3. DOS for weak coupling regime and three values 
of gate potential V. e f = 0, F/U = 0.05 and t" /U = 0.02. 
Insets show detail of Kondo peak. 

The analysis of the density of states (DOS) for various 
values of the gate potential helps to clarify the different 
circumstances involved. In Fig. 3a the DOS shows a rel- 
atively large and broad resonant peak above the Fermi 
energy which is splitted by the inter-dot interaction. It is 
the DOS of a one-body problem since the dots are empty, 
N <~ 0. The case where N ~ 1 and the conductance has 
its maximum value e 2 /h is presented in Fig. 3b. There is 
an asymmetric Kondo peak much thinner than the two 
resonant levels which are separated by the Coulomb in- 
teraction U. The inset in the figure shows the Fermi 
level located at one of the Abrikosov-Suhl sub-peaks. Fi- 
nally Fig. 3c shows the electron-hole symmetric condi- 
tion, TV ~ 2, where the Kondo resonance is just in the 
middle of the two Coulomb blockade peaks. The Fermi 
energy coincides with the position of the pseudo-gap of 
the Kondo peak, as shown in the inset. These results 
show that in the last two cases the system is in the Kondo 
regime and explain the reduction of the conductance for 
the situation corresponding to Fig. 3c relative to the one 
of Fig. 3b. 

In summary, we have obtained a numerically exact 
T = solution for the transport properties of a coupled 
two-dot system. We show that the interplay between the 
interactions that give rise to the Kondo effect and the 



anti-ferromagnetic coupling are explicitly reflected on the 
conductance of the system. The different regimes have 
been characterized through a detailed study of the charge 
at the dots and the spin-spin correlations. 
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